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Abstract: This is a design study into the capabilities of the Australian Square Kilometre 
Array Pathfinder in performing a full-sky low redshift neutral hydrogen survey, termed 
WALLABY, and the potential cosmological constraints one can attain from measurement 
of the galaxy power spectrum. We find that the full sky survey will likely attain 6 x 10^ 
redshifts which, when combined with expected Planck CMB data, will constrain the Dark 
Energy equation of state to 20%, representing a coming of age for radio observations in 
creating cosmological constraints. Keywords: methods: numerical — telescopes — galaxies: 
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1 Introduction 

With the advent of large cosmological volume galaxy 
surveys, comprised of well measured positional infor- 
mation from homogenous datasets, the measurement 
of the galaxy matter power spectrum has become al- 
most routine. The use of such power spectra in the de- 
termination of the cosmological model has been based 
almost exclusively, however, on optical techniques, e.g. 
the 2dF Galaxy Redshift Survey (2dFGRd) and the 
Sloan Digital Sky Survey (SDS^. The overall shape 
of the power spectrum is sensitive to the cosmologi- 
cal parameter constraints V = fim h and /b = fib/f2m, 
where 0.-^ and f^b are the total matter and baryon 
densities defined relative to critical, and 



h = _ffo/(100 kms ^ Mpc ^), as well as the spectral 
index of t he density fluctuations, Wg, and neutrino de n- 



sities (e.g. 


Percival et al.ll200ll:lTeemark et al. 


l2004allbl: 


ICole et al 


l2005l: lAbdalla and RawliuEsI 12007 


). Addi- 



tional information from the power spectrum can be 
gleaned by measuring the physical s cale of the so-called 
Bary o nic Acoustic O s cillat ions fe.g. lBlake and Glazebroo^ 
l2003l : iPercival et all l2010l ) either through using the 
scale as a 'standard ruler' or by combining with mea- 
surement of the Cosmic Microwave Background to break 
degeneracies. These measurements allow the determi- 
nation of the nature of the Dark Energy, through con- 
straining the equation of state parameter, w, which 
for the cosmological constant is —1. With enough suf- 
ficiently distant galaxies the variation of the equation 
of state parameter with redshift can be constrained by 
comparing this acoustic scale in different epoch s, as 
considered in, e.g. Abda lla. Blake and R awlinga (.2010l l. 

Recent advances of the speed at which radio tele- 
scopes can survey the sky to a given fiux limit point 
to the possibility of radio joining optical surveys to 
measure the matter power spectrum. The distribution 



^2dF homepage: www.aao.gov.au/2dF 
^SDSS homepage: www.sdss.org 



of these sources along the line of sight is accurately 
determined by using the redshifted emission line at 
« 21 cm of the hyperfine splitting transition in neutral 
hydrogen (HI). Previo usly surveys have been limited 
to ~ 10^ galaxies (e.g. IZwaan et al]l2005l : iLang et al.l 
120031 ') whilst the very latest HI catalogue from the 
Arecibo legacy survey, ALFALFA, is expected to find 
~ 10* objects (|Giovanelli et al.ll2005l '). In the near fu- 
ture the Chines e-built Fiv e-hundred Aperture Spher- 
ical Telescope (|Nanl |2006| ) c ould detect as m any as 
~ 10^ in the current design (|Duffv et aLlbOOSD . Ulti- 
mately however, the future for radio galaxy surveys is 
the Square Kil ometre Array, SKA0 , which may detect 
~ 10^ galaxies (|Abdalla and Rawli ngs 200^. The ini- 
tial step towards the SKA facility is a precursor known 
as the Austrahan SKA Pathfinder or ASKAjj3. The 
pathfinder consists of a much reduced number of tele- 
scopes, but still operating with a large Field of View 
(FoV) of the sky, which therefore enables the revolu- 
tionary upgrade in survey speed. 

The low-redshift precursor surveys of the SKA will 
accurately measure the properties of galaxies at low 
redshift and how these properties change as a func- 
tion of environment, for example what is the depen- 
dence of the HI mass function on local galaxy den- 
sity? The deeper precursor surveys will measure evo- 
lutionary effects. In this work we will demonstrate 
that simple estimates of the number and distribution 
of HI detected galaxies will enable ASKAP to be the 
first radio telescope to derive cosmological parameter 
constraints, able to constrain the Dark Energy equa- 
tion of state to 20%. This is a similar capability to 
previ ous optical l y bas ed measurements such as with 
2dF (|Cole et all 120051 ) but significantly more collect- 
ing area will b e needed to rival cu rrent optical surveys 
such the 6dF (iBeutler et al.ll201ll). SPSS I I luminous 
red galaxy survev fe.g. iThomas et al.ll201ll ) and Wig- 



•^SKA homepage: www.skatelescope.org 
*ASKAP homepage: www.askap.org 
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gleZ (|Blake et al.ll2O10l l surveys much less forthcoming 
optical surveys such as the SDSS III Baryon Oscilla- 
tion Spectroscopic Survey (BOS^). However we note 
that future radio surveys with, for example, the SKA 
have the capability to exceed the volume and numbers 
of galaxies spectroscopically detected compared with 
optical surveys. 

While limiting our study to the use of the power 



extent on the sky). In other words, with the large 
baselines, 2 km, available to WALLABY, most detec- 
tions are resolved and hence one must consider this is- 
sue. The positive counterpoint to this high resolution 
is that the galaxies rarely overlap within the beam of 
the telescope and the survey is effectively never confu- 
sion limited, as discussed in Section [2.21 

For an interferometer observing a resolved struc- 



spect rum in constraining cosmology fe.g. lBlake and GlazebrdoiB e one can attempt to recover some (but not all) of 



|2003| ) we note that the spectroscopic nature of radio 
surveys enable other cosmological probes. One such 
use is in measuring redshift-space distortions, a statis- 
tical measurement of infalling galaxies near large scale 
structure, which can probe the growth of structure 
on cosmic scales and hence providing a strong test of 
the validity of G eneral Relativity as well as Dark En- 
ergy models fe.g. ISong and Percival [20091 : iBlake et al.l 

mm 

over ~Mpc scales. With a relatively high den- 
sity of galaxies probing a given volume WALLABY will 
likely equal, if not exceed, constraints on this measure- 
ment by an optical survey such as 6dF (Beutler et al., 
in prep). Instead of using the redshifts one can calcu- 
late the distance to the galaxie s themselves, using rela - 
tions such as the TuUv -Fisher (iTullv and Fisheilll977l') 



the missing flux by smoothing the data, i.e. spatially 
integrating. However, some flux will be lost as spa- 
tial smoothing is akin to removing the longer antenna 
baseline pairs and therefore results in a loss of sensi- 
tivity. Exact characterisation of the effect depends in 
detail on the distribution of neutral hydrogen in each 
galaxy and the properties of the 'source finder'. This 
is a complex issue and a detailed discussion of angular 
size distributions of detected galaxies is therefore de- 
ferred to a later paper (Duffy et al, in prep). However, 
a simple approximation is that, if there are n indepen- 
dent pixels (beams) that make up a galaxy (where we 
use the example n » 1), the S/N ratio of those pixels, 
when combined, will be improved by y/n if the object 

J r ---- — ■ ■■■ is approximated by a top hat column density and ve- 

or Fu ndamental Plane (|Faber et allll987l : lDiorgovski and Da^iab itv profile. We consider this effect in more detail 
1 19871 ') ■ to estimate the Hubble flow which can then 
be removed to leave the peculiar velocities of each 
galaxy. These velocities form a velocity field, or Bulk 
Flow, that reflects the matter distribution, allowing 
constraints of the amount of mass and the fluctuation 
of ov e rdensities in the Universe (e.g.lBurkev and Tavloil 
l2004l : lAbate et al.]|2008l : IWatkins et al.ll2009l ). 

We detail the techniques and assumptions consid- 
ered in our calculation of galaxy detections in Sec- 
tion [21 in particular the effect of telescope resolution 
and galaxy inclinations in limiting galaxy counts (Sec- 
tion [53} • Utilising these assumptions we calculate the 
expected number of galaxies that the all-sky ASKAP 
(WALLABY) survey could be expected to find in Sec- 
tion [S] By constraining the matter power spectrum we 
then estimate the suitability of the WALLABY survey 
as a cosmological probe in Sectional 



in Section 12.31 and find that there is a non-negligible 
reduction in the galaxy counts of order 20% if one con- 
siders a reduction in S/N due to this effect. 

2.1 Estimating the HI signal 

As detailed in lDuffv et al.l (|2008l ) and references therein, 
the expected thermal noise for a dual polarisation, 
ripoi — 2, single beam is given by 



(1) 



2 Method 

We ha ve utilised a , signi ficantly, updated methodol- 
ogy to iDuffv et al.l (|2008l ') which analysed the poten- 
tial galaxy surveying power of the Five hundred me- 
tre Aperture Spherical Telescope (FAST). Therefore 
the reader may wish to consult that article for a more 
in-depth discussion on the following issues, including 
a consideration of evolution in the HI mass function 
(WALLABY is a shallow survey and hence likely to be 
unaffected by evolution). However there is one signif- 
icant difference between FAST and ASKAP, namely 
that the former is a single dish and the latter an inter- 
ferometer. A difference that potentially has significant 
effects in terms of resolving out extended structure. 
As we shall see the loss of signal from this effect is an 
issue for objects at all redshifts not just those closest 
to the observer (and hence with the largest angular 



for an observing time of t and a frequency bandwidth 
Au, where k = 1380 Jym^ is the Boltzmann con- 
stant and Taya is the system temperature (assumed 
to be 50 K). The effective area, Abh, calculation has 
been modified from the previous single dish calculation 
to better reflect the interferometer nature of ASKAP. 
The individual effective area of an ASKAP dish is 
the geometric area of a 12 m diameter dish, aefi, re- 
duc ed by the aperture effi ciency, expected to be Ocff ~ 
0.8 ([Johnston et al.ll2008l l. Due to computational lim- 
itations in correlating signals from all 36 dishes in 
ASKAP, WALLABY will likely use the inner core of 30 
dishes, A^dish ~ 30, which can be combined in A^porm ~ 
A'^dish(A^dish — l)/2 permutations. The resolution of the 
inner core is limited to 30" at 21cm wavelength using 
the central 2 km baselines of ASKAP. For each pair- 
wise correlation we assume a y/2 boost to the signal-to- 
noise by averaging th e real and imagin ary signal from 
a complex correlator (|Thompsonll200'7[ ) . This leads to 
an overall effective area for ASKAP of 



V2- 



kTa, 



^cffaeff ^npoiA'dish(A^diBh - l)Aut 



(2) 



^BOSS homepage; www.sdss3.org/surveys/boss.php 



where we have averaged over the complex and real sig- 
nals. 
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Typically, the beam area increases like oc (1 + 
z)^ which, if one uniformly tiles the 2 = sky, has 
the positive result that slices at higher redshift receive 
extra exposure due to the fact that observations will 
overlap. This reduces the flux limit relevant to a par- 
tic ular redshift slice by a factor (1 + z)~^, as discussed 
bv lAbdalla and Rawling^ ([200l). This is not the case 
for ASKAP however, as the number of on sky beams 
varies as a function of frequency to ensure that there 
is an approximately fixed covering area as a function 
of redshift. Hence, the flux limit for an observation. 
Slim, for a specific signal-to-noise ratio (S/N) is given 
by 



Slim = {S/N)an 



(3) 



We relate this flux to the HI mass. Mm, of a galaxy 
at redshift z in terms of the ob served flux, S, and line 
width, AVo, bv [Robert si (|1975| ) 



Mm _ 2.35 X 10" 

Mq ^ 



1 + z 



Mpc 



kms^ 



(4) 



where dL{z) is the luminosity distance to the galaxy, 
necessitating the {l + z)~^ correction for an FRW uni- 
ve rse. In a sign i ficant departure from the methodology 
of iDuffv et al] (|2008t ) we make use of the measured 
number density of objects as a function of velocity 
wi dths and HI masses d irectly from HIPASS, presented 
in IZwaan et al.l (|2010f ). With this method we auto- 
matically include the effects of angle of inclinations of 
galaxies as well as the complex velocity-structure of 
the system. 

In Fig. [1] we show the full matrix utilised not- 
ing that the histogram widths are 0.01 dex whereas 
the colour scheme is the standard number density in 
decades of mass and velocity. We emphasise that this 
represents the very latest information pertaining to the 
frequency of HI systems as a function of mass and ve- 
locity widths and, due to the limited redshift surveyed 
by WALLABY, is an ideal basis for estimating galaxy 
number counts. 

To explicitly incorporate the matrix information 
we must recast Eq. [l] first dividing by AVo (the ratio 
is termed a peak flux 5'p<=^'^ = A/hi/ AVo, implicitly as- 
suming the galaxy spectral profile is a top-hat). This 
means that for each grid cell in Fig. [1] we can also 
divide the HI mass by the velocity width (W20 is as- 
sumed to be the full velocity width of the system) and 
read off the number density of galaxies in those cells 
which lie above this peak flux. To determine what the 
minimum observable flux is (as a function of redshift) 
we substitute the flux limit of the survey (Eq. [S} for S 
in Eq. |4] allowing us to create the following inequality 
to determine when a population of HI sources will be 
detected by ASKAP 



S 



peak 



Mm 

Me 



( AVo y 

V kms-iy 



km s~^ 

2.35 X 10^ / di,{z) 
1 + z \ Mpc 



/(s/JVK 
I Jy 



(5) 



where A^ch is the number of channels that the galaxy 
will be distributed across, calculated as rounding to 




og,o [Mpc ■'] 



10 



Figure 1: We have created a matrix of the de- 
tections from the HIPASS survey as a function of 
velocity widths and inferred HI mass. This ma- 
trix currently represents the latest understanding 
in the distribution of HI detected galaxies, with an- 
gle of inclination effects as well as rotation - mass 
relations represented. 



the nearest integer number of channel AVo/dV (where 
dV in A SKAP is 4kms~^) in a calcu lation similar to 
Eq. 8 of lAbdaUa and Rawling^ l|2005f l. 

One can estimate the number of galaxies detected 
in the survey by adding the densities of all popula- 
tions that fulfil the inequality in Eq. [5] with a peak 
flux Sf°^^{z) by computing 



N{M > Afiim,z) = AQAz 



dV 

dzdO, 



dN 



dS, 



where the sky area covered is AO. and the size of the 
redshift bin is Az, dV/dzdO. is the comoving volume 
element for the FRW universe and dN / dV dS'^"''^ = 
dN/dV{dMm/dW2o) is the comoving number density 
of galaxies per unit peak flux, taken from the ma- 
trix of Fig. [1] No evolution in the HI mass-velocity 
width space has been assumed over the limited red- 
shift surveyed in WALLABY, nor is there any appar- 
ent evol ution in the in tegrated cosmic HI density out to 
2« 0.8 (|Chandl2010h or indeed z « 2 (jProchaska et all 
I2OIOI I. We calculate the average redshift of galaxies in 



the survey from N{S'P'"''' 
appropriately over z, that is. 



^ •'lim 



, z) by integrating 



/g°°zA^(SP'="'^ > Sf^f,z)dz 
j^NiSP--"^ > S^'',z)dz 



(7) 



2.2 Confusion of galaxies 

A limiting factor in galaxy surveys is the issue of con- 
fusion, whereby detections in HI are unable to be un- 
ambiguously assigned to a single galaxy. Typically HI 
surveys have previously had far greater discrimination 
between objects along the line of sight than in the 
plane of the sky. ASKAP will differ in this regard by 
enabling wide-field surveys of the sky with at least 30" 
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Parameter 



WALLABY 



Aes (m^) 

Tsys (K) 

Maximum Baseline (km) 
Angular Resolution (z=0) 

Sky Coverage (deg^) 
Total Survey Time (hrs) 

Redshift range 
Total number of galaxies 
Mean redshift of sample 
V^s atk = 0.065 /i-^Mpc 
(Tp/P at A; = 0.065 h'^ Mpc 



2668 
50 
2 
30" 
30000 
9600 
- 0.26 
673321 (848566) 
0.0492 (0.0557) 
6.9 X lO^Mpc^ 
9% 



Table 1: We summarise here the 
survey specific values of WAL- 
LABY (|Koribalski and Stavelev-SmithI l2009t) 
in addition to the strawman values of 

l2008h . 



ASKAP ( Johnston et al.l 



We consider 
the reduced baseline model for WALLABY which 
utilises the inner 30 dishes across a maximum 
2 km baseline rather than the full 36 dish, 6 km 
extent of ASKAP. We also have two numbers 
for the predicted galaxy counts, and their mean 
redshift, reflecting the effects of including the 
reduction of signal-to-noise by spatially resolved 
galaxies, as demonstrated in Fig. [2] The brackets 
ignore this effect and hence have a larger galaxy 
count. Wc consider the conservative estimate 
for ASKAP when noting the effective volume, as 
given in Equation [Ml for the typical fc-mode of 
interest (roughly the position of the first Baryonic 
Acoustic Oscillation peak). We also give the 
estimated percentage measurement error of the 
power spectrum at this scale. The maximum 
redshift we can probe the first peak out to, with 
nP = 3, is z = 0.116 with a number density of 
sources estimated at n = 1.07 x 10~''Mpc~^. 
At the second peak k = 0.125 Mpc we find 
the effective volume to be 3.6 x 10^ Mpc^ and 
ap/P = 7%. 



resolution (ASKAP baselines of 2 km) together with 
highly competitive 4kms~^ velocity resolution. It is 
therefore unlikely that confusion will play a significant 
role in limiting the number of galaxy detections in the 
WALLABY survey, an expectation we verify by mak- 
ing use of a simple analytic e stimate of the expected 
level of confusion (as given in IStavelev-Sm ith 2008). 

To estimate the confusion rates one must have a 
measure of how many galaxies are in the survey vol- 
ume, which we estimate from th e HI mass function as 
measured bv lZwaan et al.l (|2005l ). 



(p{Miii)dMm 



\ M* J 



V J V M* J 



(8) 

where a = -1.37, 6* = 1.42 x 10"2( ft"^ Mpc)"^ and 



M* = 10 Mq. The number density of galaxies, no, 
with a mass larger than Mii,n is given by the integral 
of this Schechter function resulting in the well known 
r function 



no(M > Afii„ 



eT(:..^) 



(9) 



The total HI mass contained in systems above Miim is 
given by the mass weighted integral of the Schechter 
function 



Mjii(t>{Mm)dMjii = Af*e*r 



(10) 

We are interested in understanding the confusion rates 
with sources that would significantly contribute to the 
mass of the main detection. Using Eqn.[TU]we can eval- 
uate the HI mass fraction contained in sources above 
a limiting mass; only 21% of the HI in the Universe is 
contained in systems of mass greater than M* , while 
75% of the mass is contained within systems more mas- 
sive than O.IM*. To a good approximation therefore 
we can calculate confusion rates from sources of this 
mass and greater. The typical density of sources above 
O. IM* is found using El which for the HI mass function 
of IZwaan et~all (12009 )18 0.017 Mpc"^. 

We then calculate the typical distribution of the 
galaxies on the sky, as given by the galaxy- galaxy dis- 
tribution which relates the typical density of sources 
no as a function of comoving distance r from a given 
galaxy, approximated in the non-linear regime by 



p(r) 



no 



(11) 



where ro is the correlation length. The average number 
of objects in a cylinder of comoving line-of-sight depth, 
j3, and transverse comoving radius, k, is 



n(/3,K) = 



2niip{r)dKdp . 



(12) 



For r'^ — + ffi the solution is 



n(/3, k) = tt/Sk no I 1 



7 - 



2F1 



2 2' 
(13) 

where 2^1 is a hypergeometric function. Assuming, 
reasonably, that the distribution of galaxies does not 
evolve over the redshift range probed by WALLABY 
we can make use of the HIPASS galaxy-galaxy cor - 
relation function measurements of iMever et al.l (120071 1 
(7 = —1.5 and ro = 4.7 Mpc). The cylinder diameter, 
2k, is set by the ASKAP beam. For th e natural Gaus- 
sian a ntenna distributio ns described i n lStavelev-SmithI 
(|2006l 'l and modelled in iGupta et a] (|2008l 'l the FuU- 
Width Half Maximum beam extent for WALLABY is 



= 1.4A/2000m, 



(14) 



where we have conservatively assumed that the maxi- 
mum ASKAP baseline is 2 km. The cylinder depth is 
set by the accuracy of the available redshifts for the 
confusing population, Az. 



www.pu blish . csiro. au /journals/pasa 



5 



Although the 4kms"^ velocity width of ASKAP 
will ensure that the spectroscopic redshifts of WAL- 
LABY will be measured to approximately 10~^ the 
best spectroscopic redshift estimate in practice will be 
limited by the typical Doppler width of the galaxy. 
If we assume, conservatively, that two M* galaxies 
of typical velocity width 300kms~^ are just overlap- 
ping then Az = 600kms"Vc = 0.002. Even with 
this conservative calculation the chances that there 
will be more than one galaxy in the beam volume, i.e. 
the confusion rate, is at a negligible sub-percent level 
for the mean redshift of the survey z « 0.05 (deter- 
mined in Section [3]) where the majority of detections 
lie. However, as one surveys deeper in the Universe the 
confusion rate will steadily increase as the telescope 
beam encompasses greater regions of space, further 
compounded by the lengthening observed wavelength, 
yet even at the survey edge of WALLABY, z = 0.26 
we find that less than 5% of galaxies will be confused. 

Using this formalism we can also calculate the typi- 
cal success rates in assigning an optical counterpart to 
the HI detected galaxies, assuming that such an op- 
tical catalogue contains all these HI galaxies we can 
then utilise the same number density as before. As a 
worst case scenario we further assume that only op- 
tical photometric redshifts are available , with a 'typ i- 
cal' redshift error Az « 0.05 (fHildcb ran dt et al.|[200a '). 
Thus the chances that we can unambiguously assign 
an optical counterpart (i.e. there is only one galaxy in 
the volume) occurs in more than 97% of cases at the 
average redshift of the WALLABY survey z « 0.05 
(determined in Section [S]). At the survey edge the 
success rate is 80% which is more than acceptable for 
most science cases. However, this value is a conserva- 
tive case as we could blindly assign a counterpart from 
the candidates, i.e. if there were two galaxies in the 
beam volume then we would be right 50% of the time. 
This could be further improved by using prior knowl- 
edge and choosing the largest stellar counterpart, for 
example. 

In conclusion, provided ASKAP baselines of 2 km 
are available (i.e. the most conservative case) the over- 
all galaxy number counts will be largely unaffected by 
confusion and this effect is henceforth ignored in the 
following discussion. To enable unambiguous optical 
counterparts detections for the majority of ASKAP 
HI detections we find that photometric errors of order 
Az ~ 0.05 are sufficient unless surveys probe deeper 
in redshift (by z « 0.4 the success rate is less than 
50%) or precise redshifts are needed, for example if at- 
tempting spectral stacking experiments, in which case 
a spectroscopic follow up in the optical is demanded. 



are present in the HI mass-velocity matrix. 

We estimate the angle of inclination, 9, of the galaxy 
by relating the measured velocity width AVo from the 
Hl-mass- velocity matrix, to the intrinsic linewidth width, 
AVe- The difference between the true velocity and ob- 
served velocities is proportional to sin(&). We can use 
the HI mass from the HI mass- velocity matrix to deter- 
mine the intrinsic linewidth of a galaxy, correcte d for 

broadening, using the e mpirical relation found by (|Briggs and Raol 
119931 : iLang et al.li2003l ') to be 



420 kms- 



Mh 



101" M(7 



(15) 



although we note that this relation shows a large dis- 
persion, especially for dwarf galaxies. The linewidth of 
a galaxy, AVe, which subtends an angle 9 between its 
spin axis and the line-of-sight can be computed using 
the T uUy-Fouque rotation scheme (|Tullv and Fougug 

ligssl ) 



(AKsin(e))' = {AVgf + {AVtf- 



2AVeAVt 1 - e 



2(ak; 



"e V 



6) 



AK = 120kms^i represents an intermediate transi- 
tion between the small galaxies with Gaussian HI pro- 
files in which the velocity contributions add quadrat- 
ically and giant galaxies with a 'boxy' profile repro- 
duced by the linear addition of the velocity terms. 
AVt « 20kms~i is the velocity width due to ran- 
dom motions in the disk (|Rhee and vaii Albadal[l996l : 
IVerheiien and Sancisill200ll '). 

With this definition of 9, zero corresponds to face- 
on and 9 = TT /2 to edge-on. In cases where AVg » 
AVc, one can see that AVe = AK + AVe sin 6*. For 
6' = 0, one finds that AVe = AK, in other words 
the HI dispersion in the disk, whereas for 9 = n /2 we 
recover AVe = AK + AK as expected. 

In addition there is a broadening effect, AKnst, of 
the HI profile due to the frequency resolution of the in- 
strument, R. For a range of galaxy pro files, this broad- 
ening is found to be AKnst ~ 0.557? (jBottinelli et al.l 
|199G| ). As befits a next generation ratio instrument 
the ASKAP velocity width is extremely fine, AKnst ~ 
4 kms~^, which is an insignificant source of error in 
the present discussion. 

However, for comp leteness we add A Knst linearly 
to AKfl, as argued by iLang et al.l (120031 ). to give the 
effective observed linewidth. 



AVo{9) = AK + AKn 



(17) 



2.3 Resolving rotating galaxies 

An important consideration for interferometers is the 
issue of resolving out galaxies that are larger in extent 
than the beamsize. For ASKAP, with a 2 km base- 
line this will certainly be an issue for extended sources. 
Furthermore with 4kms~i velocity resolution ASKAP 
will also probe the velocity structure of the galaxies 
themselves hence we try to estimate the size, inclina- 
tion and typical rotation velocities of the galaxies that 



The inclination angle 9 can be solved for in Eqn [16] 
by substituting the above effective observed linewidth 
('Ean ll7|l and the intrinsic linewidth from Ean ll5l We 
note that calculating the angle of inclination from the 
matrix presented in Fig. [T] in this way finds less galax- 
ies than by randomly assigning an angle of inclina- 
tion, uniform in cosine, to the galaxies (which typi- 
cally raises the completeness for WALLABY to 90%), 
in this case the more detailed calculation is the more 
conservative estimate. 
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To determine the extent of an object on the sky, 
and hence ultimately the number of beams that resolve 
the structure, we make use of an empirically derived 
relation between the HI mass of a galaxy and the HI di- 
ameter, Dhi (defined to be the region inside which the „ 

HI sur f ace density is greater than IM q pc~'^). From lBroeils and iSieel 
(llQQTl V rVerheiien and Sancisil (|200ll ') we have 



Dm ^ f Mm 
kpc ^ \10<^ »Mq 



(18) 



The o n-sky area of the galaxy can then be estimated (jMever 
I2OO8I ') using Tv{Dm/2f{B/Af where A and B are the 
major and minor axes respectively, the ratio of which 
(B/A) is equal to cos{9) which is calculated above. In 
practice we limit the smallest measurable ang le of incli- 
nation for spirals to VO.12 in accordance with lMasters et aL 
(|2003l ). We compare the apparent area of the galaxy 
on the sky, scaling by the square of the angular diam- 
eter distance dA{z), with the assumed Gaussian beam 
of ASKAP, Abeam, given by 




0.00 



0.05 



0.10 0.15 0.20 0.25 



— TvQv 



1/(4 In 2). 



(19) 



where JIfwhm was defined in Eq. [T?] previously. 

As described at the start of Section [5] we assume 
that the Signal-to-Noise of the galaxy is reduced when 
we are forced to recombine the multiple beams by 
which a galaxy is, potentially, resolved. The loss of sig- 
nal is approximated by the square root of the number 
of beams needed to cover a given galaxy, given by the 
galaxy area, Agai, divided by the beam area, Abeam- 
In practice the beam and galaxy widths are convolved 
when estimating this reduction, ensuring that even if 
the galaxy matches the beam size an additional factor 
of unity is added to this number of beams, giving a y/2 
reduction in signal in this case. Hence we reduce the 
peak flux 3'^°'^^ in the detection Eq.[5]by this geomet- 
ric factor y^l + Agai /Abeam- For the current design 
of WALLABY, with a 2 km baseline, nearly 80% of all 
galaxies are recovered after this reduction in peak flux- 

In Fig-[2]we compare the predicted number counts 
as a function of redshift for the full sky WALLABY 
survey both with (black, solid curve) and without (red, 
dot curve) the effects of resolving the galaxies. Clearly 
this effect is particularly an issue for the faint distant 
sources which are both face-on and massive to be re- 
solved out. 



3 Galaxy survey 



Figure 2: In this figure wc show the expected num- 
ber counts of galaxies in redshift bins of width 
Az = 0.01, both with and without a loss of signal 
due to the resolving out of galaxies by ASKAP 
(black solid and red dotted curves respectively). 
The effect is important, especially for systems at 
high redshift when the reduction in Signal-to-Noise 
is sufficient to push them below the detectabil- 
ity limit of ASKAP. The survey has values as de- 
scribed in Table 1. WALLABY is approximately 
80% complete for a baseline of 2 km. On the right 
axis, in blue, we plot the limiting HI mass as the 
dotted blue curve for a signal to noise detection 
of 5ct in one pointing in redshift bins of width 
Az = 0.01 and a velocity width of 200kms~^. 



The expected number counts as a function of red- 
shift on completion of the proposed survey is shown 
in Fig. [2] as the solid black curve with the actual total 
number of detections and mean redshift of WALLABY 
given in Table 1. 

4 Cosmological Parameters 

Using the predicted galaxy number counts for WAL- 
LABY we can estimate the errors on the galaxy power 
spectrum at the mean redshift of the survey z — (z) ^ 
0.055 and ultimately the expected constrains on cos- 
mological parameters. P{k, z) is related to the power 
spectrum P{k, 0) by 



Diz) 



In this section we combine our estimates of the de- 
tectability of galaxies from the previous section^ with 

the ASKAP strawman flgures (|Johnston et a l.' 2008) 

and the specifics of the WALLABY survev (jKor ibalski and Stavelev-SmithI ^ 
l2009l l. as summarised in Table 1. In Fig. [2] the dashed n 
blue curve indicates the expected neutral hydrogen 
mass limit as a function of redshift in redshift bins 
of width Az = 0.01 for a single pointing of ASKAP. 
The redshift depth of WALLABY is such that the sur- 
vey ends when the mass limit approaches < 10^^ M0, 
which is the apparent maximal limit of HI systems. 



P{k,z) = [bD{z)fP{k)., 



(20) 



where b is the bias parameter and D{z) is the growth 
factor computet^ from 



-E{z] 



{l + z')dz' 
[E{z')Y^ 



(21) 



where E{z) = H{z)/Ho and by construction D{z = 
0) = 1. 



^ Using the excellent publicly available icosmo pack- 
age l IRefregier et al.ll2011l '). 
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Errors on the power spectrum are due to two fac- 
tors: sample variance, i.e. the fact that not all k modes 
are measured, and shot-noise, which is the effective 
noise on the measurement of an individual mode. The 
total error ap on the measurement of the power spec- 
trum, P(k,z), fo r a given fc with bin width Afc can 
be expressed as dPeldman et al.1 1 19941 : iTegmarkI Il997l : 
iBlake et al.ll2006D 

where P = P(t.z); „ = n{,) = j-^^^^^^dM 

is the number density of galaxies which are detected 
(making nP dimensionless) and m is the number of 
k-modes in a survey of total volume V , with m = 
2Tvk^Ak V/{2tt)^. The ability of a survey to probe cos- 
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Figure 3: The top panel shows the power spec- 
trum (black curve) for the WMAP7 cosmol- 
ogy ( Komatsu et al.|[2009() . at the mean redshift 
of the WALLABY survey. The points are the ex- 
pected errors on the measurement from the Stt sky 
survey with WALLABY. The bottom panel is the 
same curve, normalised to a reference run with no 
baryons. The vertical dotted line indicates the 
smallest fc-value (fc = 0.15/i^^Mpc) we consider 
for analysis. 



Surveys 


w 




Ho 


Planck 


-1.02 ± 


0.28 


73.3 ± 10.4 


Planck + WALLABY 


-0.98 ± 


0.12 


71.2 ±4.2 


Planck + 2dF 


-1.00 ± 


0.13 


72.0 ±4.6 



Table 2: Shown are the predicted cosmological 
parameter estimates when projected Planck CMB 
data is used alone, then in combination with WAL- 
LABY and also when used with an existing optical 
survey 2dF. Including the matter power spectrum 
measurements results in a factor two improvement 
in the constraints on the Dark Energy equation 
of state parameter w and the Hubble constant 
Hq. Note that we performed a 6 parameter cos- 
mological fit but that Planck has such small er- 
rors on most parameters that a survey with less 
than a few 10^ sources is unlikely to improve the 
estimates, with the exception of w and Hq, and 
for completeness we list the other variables which 
were unchanged when 2dF or WALLABY sur- 
veys were included; [flhh'^ ,^ch^ , rig, \og{lQ^° As) ,t] 
were constrained to be [0.0227 ± 0.0002,0.1099 ± 
0.0015,0.964 ± 0.005,3.06 ± 0.01,0.092 ± 0.006]. 



mological parameters can be estimated, and compared, 
with the effective survey volume VcB given by 



Voff (fc) = An 



dV 
dzdQ 



{z)dz . (23) 



As argued by ISeo and EisensteinI 
great gain in probing beyond nP 



(120031 1 there is no 
= 3 hence we limit 
the volume of our survey to a maximum redshift Zmax 
when n{z — Znia.yi)P{k) = 3; conservatively assuming 
a constant number density within this volume, hence 
the effective volume is now evaluated as 



vMk)^An 1 + 



1 



n(Zmax)-P(fc) 



dV 
dzdQ. 



{z)dz , (24) 



evaluated in Table 1 for two fc-modes of interest (fc = 
0.065 and 0.125 /iMpc~^), approximately correspond- 
ing to the baryonic peaks in the power spectrum. Note 
that we have assumed a constant weighting function 
for the galaxies used in the power spectrum measure- 
ment, a so-called 'number-weighting' scheme. With 
our assumptions of constant number density within a 
given maximum redshift the WALLABY survey is es- 
sentially a uniform survey (i.e. with a window function 
which is effectively the identity matrix) and the covari- 
ance matrix of the k-bands can be well represented by 
a diagonal covariance mat rix with elements equal to 
(4/3)^(PV^) as argued in lBlake et al] (120061). 



In this work we create a power spectrum (|Lewis et al.l 
I2OOOI ') based on t he latest WMAP7 Ma ximum Likeli- 
hood cosmology (|Komatsu et al.l 120091 ) for the stan- 
dard ACDM model with values [Q^, Qui, Qa, h, w, ris, as] 
given by [0.0451, 0.271, 0.729, 0.703, -1, 0.966, 0.809]. In 
the top panel of Fig. [3] we demonstrate the power spec- 
trum with the expected errors from WALLABY. In 
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the bottom panel of this figure we have normahsed the 
matter power sp ectrum by a reference n o-oscillation 
power spectrum (jEisenstein and Hulll998| ^ to aid visu- 
alisation of the peaks. It appears that with the er- 
rors from WALLABY it will be challenging to iden- 
tify the 'baryonic wiggles'. To quantify this we have 
performed a test on the no-oscillation model with 
the datapoints in Fig.|3]to determine if the datapoints 
have sufficiently small errorbars to rule out this wiggle- 
less model. The x^/d-o.f value is 2.2/7, this means 
that WALLABY will only detect the BAO peaks with 
\/2.2a significance. Although not conclusive this anal- 
ysis does suggest that WALLABY will not significantly 
detect the Baryonic Acoustic Oscillations (BAO), there- 
fore we do not attempt to use these oscillations as stan- 
dard rulerfl. 

In accordanc e with our previo us method of the 
analysis of FAST (|Duffv et al.ll2008l l we limit our inves- 
tigation of the power spectrum to band- powers over 
the range 0.005 < k/hMpc~^ < 0.15. The maximum 
wavenumber chosen is a conservative cut on the power 
spectrum to ensure that we pe rform our analys is in the 
linear regime (as suggested in lCole et al.ll2005h . It has 
been suggested (Beutler et al, in prep) that an HI sur- 
vey such as WALLABY could safely trace even smaller 
scales (~ 0.2 h Mpc"^). This is because of the inherent 
low bias of HI detections which are tidally stripped in 
high-density regions, hence a blind HI survey will nat- 
urally avoid these regions, ensuring that the nonlinear 
effects due to peculiar velocities of galaxies in these 
groups and clusters will be lessened, effectively extend- 
ing the linear regime to smaller scales. For the Cosmic 
Microwave Background (CMB) data we have assumed 
that we have full polarisation information for Planck 
by considering the temperature T and E-type polarisa- 
tion anisotropies for I < 2400 (including cross spectra) , 
and assumed that they are statistically isotropic and 
Gaussian. The noise in the CMB data is also assumed 
to be isotropic and is based on a simplified model with 
Ntt = Nee /'I = 2 x 10~"'^K^, having a Gaus sian 
beam of 7 arcminutes (jPlanck Collaboration|[200^ '). In 
the Markov-Chain Monte Carlo analysis we analyti- 
cally marginalise over the bias parameter b in Eq. [20] 
and assume that the priors around each cosmological 
parameter are flat, with a width safely outside that 
aUowed by WMAP (|Komatsu et al.ll2009l '). 

The values shown in Table 2 are best-fit cosmo- 
logical values (|Lewis et al.ll200"^ for a variety of dif- 
ferent parameters using expected Planck CMB data 
alone and then combined with WALLABY (or 2dF as 
a comparison) . By combining CMB data with the mat- 
ter power spectrum measurement, which itself isn't a 
function of w, we break the degeneracy between w and 
h that occurs when calculating the distance to the sur- 
face of last scattering from the CMB. Thus the main 
effect of ASKAP is to reduce the error on h and w by 
a factor two on the value achieved with Planck alone, 
as demonstrated in the top panel of Fig. [4] In the 
bottom panel of Fig. [4] we compare WALLABY with 

'^In lBeutler et al.l (120111) they also predict that the BAO 
peak is marginally detectable in the correlation function 
using WALLABY, at 2.1 ±0.7(T which is in agreement with 
our estimate 
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Figure 4: Cosmological parameter error ellipses for 
degeneracies between the Hubble constant and the 
equation of state of Dark Energy in the top panel 
and the matter density in the Universe F = fim h 
versus the universal baryon fraction /b in the bot- 
tom panel. In the top panel we consider the case 
of a CMB measurement alone from Planck in blue 
and in yellow the improved constraint from the 
measurement of the power spectrum by WAL- 
LABY (2dF gives a very similar result). In the 
bottom panel we compare the results from the 2dF 
spectroscopic survey in red and the WALLABY 
estimates in yellow again. With these expected 
results ASKAP will likely be the first radio tele- 
scope to successfully measure cosmological param- 
eters from the matter power spectrum. 



an existing optical based measurement of the matter 
power spectrum from 2dF, the two error ellipses are 
of similar size, graphically illustrating the parameter 
estimates in Table 2 that WALLABY will be the first 
radio telescope to infer cosmological parameters from 
the matter power spectrum. 
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5 Conclusion 



Blake, C. and Glazebrook, K. 2003, Ap. J. 594, 665 



As is clear from the galaxy survey estimates for WAL- 
LABY, ASKAP will likely be the first radio telescope 
to measure the matter power spectrum, constraining 
cosmological parameters to the level attained by the 
optical survey 2dF in 2005, this will represent a com- 
ing of age for radio astronomy. To match, and ulti- 
mately surpass, current surveys such as WiggleZ will 
demand the full version of the SKA. In creating a radio 
based galaxy survey the matter power spectrum will be 
analysed using a different galaxy tracer together with 
different survey selection effects than the 2dF spec- 
troscopic optical catalogue. Although not considered 
here there is the additional possibility of using the ve- 
locity field of the galaxies to gain additional cosmolog- 
ical constraints. As regards to a full local sample of 
~ 6 X 10^ HI detected galaxies the science case is in- 
triguing for the determination of star formation in the 
local Universe. When coupled with deeper surveys on 
ASKAP, such as DINGO, that can determine the evo- 
lution in redshift out to z = 0.4, the combined output 
will be a significant dataset for years to come. 
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